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Abstract: 

Purpose — The purpose of this paper is to develop a serial redundant manipulator system applied in nuclear fusion 
environment. It will allow remote inspection and maintenance of plasma facing components in the vacuum vessel of 
fusion device without breaking down the ultra-high vacuum condition during physical experiments. 
Design/methodology/approach — Firstly, considering the dynamic sealing of actuators to avoid polluting the vacuum 
condition inside fusion reactor, the mechanical design of robot system has been introduced. The redundant manipulator 
system has 11 degree of freedoms in total with an identical modular design. Besides, to improve the position accuracy, 
an error prediction model has been built based on the experimental study and back-propagation neural network (BPNN) 
algorithm. 

Findings — Currently, the implementation of the manipulator system has been successfully finished in both atmosphere 
and vacuum condition. The validation of BPNN model shown an acceptable prediction accuracy (94%~98%) compared 
with the real measurement. 

Originality/value — This is a special robot system which is practically used in a nuclear fusion device in China. Its design, 
mechanism and error prediction strategy have great reference values to the similar robots in vacuum and temperature 


applications. 
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1. Introduction 


Experimental Advanced Superconducting Tokamak 
(EAST) is the world's first fully superconducting tokamak 
fusion device with non-circular cross-section which was 
built in China (Wan et al., 2006). In recent years, with the 
increasing device performance and_ experimental 
parameters, EAST tokamak has made a series of important 
research results and scientific discoveries (Li et al., 2013; 
Wang ef al., 2014). However, EAST inner components 
will also be facing increasingly tough operating 
environment with huge heat flux as the plasma current 
during discharges becoming higher and higher. It leads the 
plasma facing components (PFCs) easily to be damaged, 
in which case the effective running time cannot be 
guaranteed (Yang et al., 2016; Li et al., 2016). Therefore, 
it is essential for the timely maintenance based on the 
condition of damaged internal components in the 
experimental period. 

The AIA robot is a related application which was 
specially developed for ITER fusion reactor (Shimomura 
et al., 1999) and demonstrated in Tore-supra tokamak 
device (Gargiulo et al., 2009). It is able to inspect the in- 


vessel components under 120°C and 10° Pa condition 


with a 200 mm maximum position error due to the 
flexibility of the 10-meter cantilever structure (Perrot et 
al., 2012). The poor position accuracy means that the 
robot can only run in the middle of the vacuum vessel 
rather than close-range operate or even touch the tokamak 
wall from a security point of view. 

For the purpose of real-time detection and rapid 
repairing of the damaged internal components during 
plasma discharges, EAST articulated maintenance arm 
(EAMA) system has been developed since 2013. It will 
allow remote inspection and simple repair of plasma 
facing components (PFCs) in EAST vacuum vessel (VV) 
without breaking down the ultra-high vacuum condition 
(~10°Pa). Due to its long-reach mechanisms with a 
weight more than 100 kg, similar with the AIA robot, the 
gravity effect will cause huge flexible deformation, which 
is unacceptable for running inside a narrow and complex- 
shaped space as EAST VV. Figure | shows the cantilever 
structure of fully assembled EAMA robot without end- 
effector and the finite element simulation results for the 


position errors. 
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Figure 1 Cantilever structure of EAMA robot 
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Note: (a) The fully assembled EAMA robot without end-effector; (b) The finite element simulation results for the position errors of 


different links. 


Elastic deformations due to gravity will cause robot 
trajectory planning to be quite a complicated issue, 
especially for serial long-reach manipulators. Error 
prediction strategy should always be studied in advance if 
high accuracy was required. In fact, research on the 
robotic flexibility has become a matter of great concern in 
recent years because of the pursuit of high accuracy and 
reliability. Normally, two kinds of methods could be 
considered to build the flexible model for error predicting. 
On one hand, the mathematic model of flexible multibody 
dynamics can be derived from classical mechanics 
theories such as Newton-Euler formulation, Lagrange 
formulation, etc. Then either the finite element method 
(FEM) or assumed modes method (AMM) can be used to 
truncate the dynamic formulations for rapid solution with 
an appropriate accuracy (Dwivedy and Eberhard, 2006). 
On the other hand, for some complex systems with several 
uncertain coupling factors, the flexible model can be 
identified using some computing algorithm without 
considering the complicated mechanical model and highly 
nonlinear formulations, such as fuzzy algorithm (Schiavo 
and Luciano, 2001), genetic algorithm (Chang, 2007), 
artificial neural network approach (Liu, 2012) , etc. 

For the case study of EAMA manipulator, the 


mathematical modelling using classical mechanics 


theories is difficult to obtain high accuracy and efficiency 
owing to two reasons. Firstly, the complicated mechanism, 
transmission chain and multi links make the accurately 
flexible modelling to be quite a difficult issue. Secondly, 
the dynamic formulations with high nonlinearity is not 
necessary as the manipulator speed is extremely low(less 
than 0.5 degree per second), which means that the 
dynamic behaviors are slight enough compared with the 
gravity effect. Therefore, static modelling can be a better 
choice for error predicting with a relatively high accuracy. 

In this paper, the conceptual design of the manipulator 
system was firstly introduced. Then, an experimental 
study was deployed to measure the errors of EAMA 
prototype assembled with different loads through a load- 
deflection platform. Based on the experimental data, a 
static error prediction model was built and trained by 
utilizing back-propagation neural network (BPNN) 
algorithm. The results show an acceptable prediction 
accuracy while around 5% predicting error compared with 
the real measurement. Furthermore, the mathematic 
formulation for static loads of robot joints in arbitrary 
position were derived. The calculated joint loads can be 
treated as the input of the trained BPNN model to predict 
the robot errors in arbitrary positions and postures without 


any experiments and measurements in future. 
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2. Manipulator system design 


EAMA system consists of a highly redundant manipulator: 


one mobile base, 5 serial arms and an end-effector (CCD 
cameras, gripper, etc.) for dedicated functional operations. 
Besides, a storage cask facility has been developed to 
maintain the equivalent vacuum environment for the 
manipulator before docking with EAST vacuum vessel. 
Figure 2 shows an overall schematic view of EAMA 
system and the robot specifications is given in Table 1 (Shi 
et al., 2016). 


Figure 2 Layout of EAMA system 
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Table 1 Specifications of EAMA manipulator 


Vacuum 10° Pa to 1 atm 
Temperature Running:80 °C; baking:120 °C 
Workspace In-vessel: R=1920, a=550 mm 
DOFs 1(base) + 7(arm) + end-effector 
Payload 25kg 
Dimension Radius: 160mm, length: 8.8m 
Weight <100 kg (arm) 
Function Inspection and simple repair 
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For the manipulator arms, an identical modular design 
with a parallelogram mechanism has been adopted for all 
of the 5 arm segments. Each segment can provide the 
motions of both rotation and elevation by integrating yaw 
and pitch joints inside one modular arm segment. As 
shown in Figure 3, the parallelogram structure is 
composed of two clevis (yaw joints), horizontal rods, arm 


Figure 3 Modular design of EAMA robot arms 
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tube and diagonal rod (pitch actuator). The five-bar 
mechanism can produce huge reduction ratio and 
withstand strong torques generated from the cantilever 
arms and the gravity effects. Besides, the axis of yaw joint 
will be always vertical owing to the motion characteristics 


of the parallelogram structure. 


Titanium robot tube 
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Note: Five segments with identical principle design, the 1* to 3™ segments have only rotation motions by using rigid rods instead of 


elevation actuators; the 4" and 5" segments have both rotation and elevation motions. 
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The motion actuators are placed inside robot tubes: 
rotation actuator connected to the tube while elevation 
actuator located in diagonal of the parallelogram. To avoid 
polluting the ultra-high vacuum condition inside EAST 
vacuum vessel, the motion actuators with the 
consideration of dynamic sealing have been developed. 
All high-speed lubricated by high temperature grease are 
sealed by SS bellows while the All high-speed 
components (motors, gears, etc.) which should be 
lubricated very well, are sealed in boxes by welding 
bellows so that high temperature grease can be used. 
Meanwhile, some low-speed components (joint bearings 
and bushes) are considered to use MoS2-Ti-C coating 
films for solid lubricating. The detailed design of motion 
actuators can be summarized as follows: 

(1) yaw actuator: 

As shown in Figure 4, firstly two planetary roller screw 
divide the rotation motion produced by high-temperature 
motor and reducer into two parallel linear motions through 
a gear group. The two screws will move at a same speed 
but opposite directions due to the opposite threads of 


screw nut. Then, the screws are connected with two steel 


Figure 4 Detailed design of EAMA yaw actuator 
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Figure 5 Detailed design of EAMA pitch actuator 
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cables which can deliver the linear motions to yaw joint 
through bellows welded together with the seal box. The 
cables will be finally assembled with a pulley system 
attached to yaw joint to produce the rotation motion. The 
whole transmission process can be summarized to be a 
“rotation-linear-rotation” chain. Two benefits can be 
provide by this chain: dynamic sealing to protect the 
vacuum condition and huge reduction ratio to generate 
enough driving torque (1:30820 from motor to yaw joint). 
(2) pitch actuator: 

As shown in Figure 5, the pitch actuator has a similar 
“rotation-linear-rotation” chain with yaw actuator. The 
difference is only one roller screw is used here to transfer 
the rotation to linear motion. As the pitch actuator is 
located in the diagonal position, the changes in diagonal 
length will lead the whole parallelogram structure to do 
elevate while other links’ length are fixed. The reduction 
ratio can reach up to 1:51660. 

Currently, all of the manipulator components have been 
developed. The whole EAMA robot system was also 
successfully implemented in EAST vacuum vessel 
(Figure 6). 


SUS316L Bellow 
(Thickness 0.1 2mm, d=16-25mm 
Operation length 50-190mm) 


Seal Box 
(Thickness 0.5mm) 


Motor and Reducer 
(Faulhaber 3257G024CR; 
Faulhaber 32/35, 1:134) 


Free space 
(113mm, 105mm needed fort45°) 


Brake 
(Faulhaber MBZ 24V) 


Reducer 
(Faulhaber 32/3S, 1:246) 


Extend Rod 
(M10, Grade9.8,14.9kN) 


Roller Screw and Nut 
(GSA RGTB 15.2, Co=24kN 
d=15mm, p=2mm, stroke 135mm) 


201710.00095v1 


chinaXiv 


Figure 6 The implementation of EAMA manipulator system 
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Figure 7 The illustration and implementation of load-deflection platform 
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Note: different payloads were attached to two adjustable beams in x and z direction to generate different torques on robot joint. 


4. BPNN Prediction model 


In a BPNN algorithm, several strongly coupled neurons 
are used to approximate the nonlinear functions. The 
learning process can be summarized as two steps 
(Negnevitsky, 2005). Firstly, a set of training data is input 
to the multi-layer network to obtain relevant output. Two 
kinds of functions will be used: a linear transfer function 
with different weights and thresholds (Eq. 1) and a 
sigmoid activation function (Eq. 2). Secondly, an error is 
calculated by comparing the output with its expected 
value. Normally the expected values are obtained from 
experimental measurements. Then the error is propagated 
backwards through the network to modify the weights and 
thresholds of the coupled neurons. The rule for modifying 
weights and thresholds should always aim at reducing the 


error. 
X =} xw, -0 (1) 


1 
Lre” = 


Y sigmoid __ 


For error prediction of EAMA manipulator, a final 
three-layer network has been built as the topological 


structure shown in Figure 8: 


Figure 8 The topological structure of BPNN model 
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(1) 4 neurons in input layer: 0 — the pitch angle of robot, 
T; — the torque along z axis, T; — the torque along x 
axis, Fy — the equivalent gravity; 

(2) single hidden layer with 8 neurons; 


(3) 3 neurons in output layer: the errors of end clevis in 


different directions: Ax, Ay and A@, ; 
(4) wy and 6; denote the weights and thresholds between 
input and hidden layer; w; and 0; denote the weights 


and thresholds between hidden to output layer. 


The 272 sets of sample data were divided into two parts: 

° 243 sets of data were utilized for network training, 
among which 70% for training, 15% for testing and 
15% for validation; 

$ 29 extra measurement data were utilized to check the 
prediction accuracy of the BPNN model after 
training. 

The mathematical model was established by the Neural 
Network toolbox in Matlab environment. The network 
training function is Traingdx, which is able to 
automatically adjust learning rate based on the classical 
BP algorithm (Shi et al., 2009). The training was ended at 
epoch 6751 while the mean squared error was converged 
to 6.7717e-5 (Figure 9). 


Figure 9 The training process of BPNN model 
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To evaluate the prediction performance of the trained 
BPNN model, the extra 29 sets of load data were 
calculated by the offline BPNN model with the trained 
weights and thresholds listed in Table 2 and Table 3. 
Figure 10 gives the fitting situation between predicted and 
measured values which indicate a prediction error range 
from 2% to 6% (Figure 11). Considering the huge 
of EAMA 


manipulator, the predicting accuracy is acceptable. The 


flexibility and complicated structure 
trained BPNN model can be integrated the control system 
to compensate the deflection error and improve the final 


position accuracy. 


Figure 10 The fitting situation between predicted and Figure 11 The prediction error of the trained BPNN 
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Table 2 The trained weights and thresholds between input and hidden layer 


j wi (i=4,j=8) O; 

1 -0.22532639 -0.021457862 0.000143862 0.021266973 1.168193494 

2 -0.36576382 0.002737846 0.000151412 0.061871295 0.506266839 

3 0.374547274 0.02131777 -0.000453784 -0.022556882 0.265934803 

4 -0.324228676 0.000747825 0.000236239 0.016785358 0.113202306 

5 -0.313983181 -0.00378339 0.000254887 0.005711087 -0.02688637 

6 -1.88891379 -0.004668008 -0.0000337 -4.908188396 -2.751156441 

7 0.389308685 -0.011864381 -0.000147494 0.045723211 0.828719931 

8 0.33108 -0.0137 -0.00012 0.036749 0.816560631 

Table 3 The trained weights and thresholds between hidden and output layer 

k Vik G=8,k=3) Ok 
1 46.3462 -19.279 -17.9038 54.933 -81.8616 -0.01364 89.21414 -121.38 -11.933 
2 7.83777 -30.161 -15.0365 181.49 -183.839 -0.05781 -73.7377 71.93208 -12.3189 
3 -23.872 9.81833 -19.0517 30.596 -69.5353 0.016233 106.448 -146.954 41.26604 


5. Formulation of static loads 


Since the motion speed of EAMA manipulator is 
quite slow (less than 0.6 deg/second), the 
dynamical characteristics are not significant 
when considering the payloads on each joint. To 
simplify the calculation, static loads due to 
gravity were derived using the model built in 
Figure 12. With the load results, the deflections of 
all segments can be predicted with the trained 


BPNN model in former section. 


Figure 12 Modelling of static loads 
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Considering the static loads on the end point O; of 
Link i, firstly the gravity load can be easily 


written as Eq. 3: 


= 3) 


For the torques, the effect on Link i caused 
by Link j was firstly considered as Eq. 4, among 
which the torque was divided into two parts that 
respectively represent the torque along x (M a ) 
and y axis ( M p) in the local coordinate system 
of Link CS (i). Besides, O,P; © is the 


position vector of point P; (mass center of Link j) 
with respect to point O; in CS (i) while 


g,=[0 g 0 oj and g =[|g 0 0 o} 


7 > p’ 
= ie | =| e (4) 
’ m,g,-0,P, 


Furthermore, OF © can be derived into Eq. 5, 


where 7 is the transfer matrix from CS (i) to CS 
(j) which can be obtained using D-H method 


(Denavit and Hartenberg, 1955). r is the 


position vector of mass center P;in CS (j), which 
is determined by the geometry and mass 


properties of Link j. 
OP. =rir; 
itj Tti Fj (5) 


Combining Eq. 4 and Eq. 5, the total torque 
applied on the end point O; of Link i can be 


derived as: 


= (6) 


n 


X (ms, TiTi) 


j=i 


wepe 


With the formulations above, the static loads 
on each joint can be calculated when the 


manipulator position is known. The loads can be 


treated as the input of the trained BPNN model 
for error evaluation without any experiments and 


measurements in future. 
6. Conclusions 


An articulated manipulator system applied in 
fusion environment has been developed in China 
for the purpose of real-time detection and rapid 
repairing of the damaged internal components in 
EAST tokamak device. The mechanical design of 
EAMA system was introduced firstly in this paper, 
including the modular parallelogram mechanism, 
the “rotation-linear-rotation” actuator design for 
dynamic sealing etc. Besides, to predict and 
compensate the flexible errors due to gravity 
effect, a load-deflection platform directed to 
EAMA prototype was built. Based on the 272 sets 
of deformation data, a BPNN model was 
established by the Neural Network toolbox in 
Matlab environment. After training, the fitting 
BPNN predicted and 


measured 


situation between 


experimental values indicate a 
prediction error range from 2% to 6%. Finally, the 
static loads of each manipulator link was 
formulated, which can be integrated with the 
BPNN model for the error evaluation in control 
system. The conceptual design and error 
prediction strategy introduced in the paper can 


give beneficial reference to the similar robotic 


applications in vacuum and temperature 
condition. 
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